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Abstract 

The cosmic ray data of PAMELA/ATIC may be explained by dark matter decay with a decay 
rate T^\f ~ 10~^^ sec^^ ~ 10^^^ eV, an energy scale which could not be understood within the 
framework of the standard model or its simple super symmetric extension. We propose anomaly 
induced dark matter decay to exponentially suppress the decay rate, and apply to a supersymmetric 
extension of the Ma's inert Higgs model of the radiative seesaw mechanism for neutrino masses. In 
this model the lightest right-handed neutrino ipN and the lightest neutralino x can fill the observed 
necessary dark matter relic, and we find that tp^ can decay into x through anomaly with a right 
order of decay rate, emitting only leptons. All the emitted poistrons are right-handed. 

PACS numbers: 



I. INTRODUCTION 



Recent astrophysical observations [l| and neutrino oscillation experiments [21] require an 
extension of the standard model (SM) so as to include dark matter as well as to incorporate 
a generation mechanism for small neutrino masses. At the moment, however, we know about 
dark matter only a little; the constraint on its mass and abundance, but nothing about its 
detailed feature is known. Consequently, there are many consistent models for dark matter. 
Representative possibilities may be summarized as follows: 

(i) Dark matter is a stable thermal relic, so that its relic abundance and annihilation cross 
section are strongly related to each other. The lightest neutralino in supersvmmetric models 
with the conserved -R-parity is a well studied example of this category [3!]. Another well 
motivated example may be a stable neutral particle in the radiative seesaw scenario 0], 
which is an alternative model of the seesaw mechanism to generate small neutrino masses 
js]. In fact, there exits similar models and the nature of the DM candidates in these models 
has been studied 



(ii) Dark matter consists of multiple components If some of them are unstable, dark 
matter can contain thermal components as well as non-thermal ones which can be produced 
by the decay of unstable components. In this case the dark matter relic abundance at present 
and the annihilation cross section of the dominant component need not to be related. 
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20|. 



(iii) Dark matter is not stable and is decaying with a very long lifetime 

In any case it will be crucial for the study for going beyond the SM to know which class the 
true dark matter model belongs to. Since the above mentioned dark matter models predict 
different signals for the annihilation and/or decay of dark matter in the Galaxy, it may be 



possible to use the data from these observations to distinguish the dark matter models |21 |. 



The positron excess in the recent PAMELA observation 22| and ATIC data of e"*" + e flux 



231] are such examples. PAMELA data show a hard positron excess compared with the 
background but no antiproton excess, while ATIC data show the excess of e"*" + e~ flux at 
regions of 300-800 GeV. A lot of works have been done to explain these data within the 
framework of dark matter models (see, for instance, 26| and references therein). However, 
the observed positron flux requires much larger annihilation cross section or much larger 
dark matter density than the ones needed for the explanation of WMAP data. The required 
large factor in the latter feature is parametrized as a boost factor in the references^. So, it 
seems very difficult to give a natural explanation for the boost factor for the type (i) dark 
matter. 



In this paper, following 27|], we consider a supersymmetric extension of the radiative 
seesaw model for the neutrino mass to understand the data obtained by the PAMELA and 
ATIC experiments. The radiative seesaw model is attractive in two respects: (a) The non- 



^ The necessary enhancement for the s-wave annihilation can be partly covered by the Sommerfeld effects 



24l | or others [2 



2 



vanishing small neutrino mass and the presence of a dark matter candidate are closely related 
through a discrete symmetry Z2. (b) The dark matter candidate in this model couples only 
with leptons but not quarks. This feature is favorable for the above mentioned PAMELA and 
ATIC data. However, a large boost factor still has to be introduced to explain the observed 



positron flux in this model 28|, |29|, |30|] . In our supersymmetric extension of the model this 



problem is overcome as follows. There are two kinds of stable neutral particles corresponding 
to two discrete symmetries, R and Z2, where R is the R parity in supersymmetric theories, 
and Z2 is mentioned above. If one of these discrete symmetries is broken, the heavier one 
can decay to the lighter one. We propose that this breaking can be induced by anomaly 
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to realize an exponentially suppressed decay rate of the heavier dark matter. 
It should be noted that the smallness of this decay rate is a crucial ingredient for the 
explanation of the observed + e~ flux. Moreover, due to the very nature of the model, 
only lepton pairs can be produced through the dark matter decay. We show that both data 
of PAMELA and ATIC can be described well simultaneously in this scenario. The model 
for dark matter proposed in this paper gives a concrete realistic example of type (ii). 



II. ANOMALY INDUCED DARK MATTER DECAY 



The stability of the dark matter is usually ensured by an unbroken discrete symmetry Z. 
If the discrete symmetry is broken, the dark matter can decay. The preferable decay modes 
depend on how Z is broken. However, its life time will be too short tdm ^ (S/7r)rnivr>M — 
10~^^ sec for rriDM — 1 TeV, unless the Z breaking is extremely weak [12]- 20|. Such 
suppression may occur if Z is broken by GUT or Planck scale physics [l6|-|2(' 

Here we would like to suggest an alternative suppression mechanism which is based on the 
observation that if a symmetry, continuous or discrete, is anomalous, then non-perturbative 
effects can generally induce non-invariant terms, like quark masses in QCD. Although the 
discrete symmetry Z in question can be anomaly free with respect to the SM gauge group, it 
can be anomalous at high energy when imbedded into a larger discrete group, because heavy 
particles can contribute to discrete anomalies 
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3J] . If the discrete symmetry is anoma^ 
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at high energy, non-perturbative effects can produce e~^^^^ in the superpotential 3ll.l32l.l33|. 
which is Z invariant. This is because the dilaton superfield S transforms inhomogeneously 
under the anomalous Z, where $ is a generic chiral super field, and 6 is a certain real number 



see also [37|, l38|). Below the Planck scale, where the dilaton is assumed to be stabilized at 



a vacuum expectation value of 0(1), the factor e"'"^'^^ can work as a suppression factor for 
the noninvariant product 

Let us estimate the size of the suppression. To this end, consider a (chiral) Z^ symmetry 
in a gauge theory based on the gauge group G and assume Z^ is anomalous. Then the 
Jacobian J of the path integral measure corresponding to the Zjv[G']^ anomaly can be written 
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as [35|,y,|38| 

J = exp (-^Ag I d'x 2Aix)^ ,A{x) = ^.e^"'" F^^F^^, (1) 

where F^j, is the field strength for G. Since the Pontryagin index J d*x A{x) is an integer, 
AQ = mod N/2 means anomaly freedom of Z^- In the anomalous case, we have AQ = k/2 
with an integer k < N. (So, AQ/A^=l/4 for an anomalous Z2, for instance.) This anomaly 
can be cancelled by the Green-Schwarz mechanism 36|], which defines the transformation 
property of the dilaton supermultiplet 5* = {(p+ia, ips, Fs), where ip (a) is the dilaton (axion) 
field, and they couple to the gauge field as 

= -| f;.^'^'^^ - ^e'^"''' (2) 

To cancel the anomaly ([1]), the axion a has to transform according to a — > a—{l/2TT){AQ/N). 
Therefore, the Z^v charge of exp(— feS*) becomes C if 6 = Atc'^C/AQ. Since < if >= l/g"^ ~ 
0(1) at the Planck scale, the expression exp{—bS) would then yield a suppression factor SF 
such as 

SF ~ exp(-47r2C/AQ), 
{SF f ~ 10"^^ 10"^^ 10"^^ for C/AQ = 8/5, 2, 10/4, (3) 

where C and 2AQ are defined modulo ^. 

Do we need such a big suppression? According to HQ, to explain the PAMELA/ATIC 
data, the decaying dark matter should decay with a life time of ~ 10^^ sec, which cor- 
responds to a decay width Tndm ~ 10"^^ x (ITeV) ~ 10"^° x (lTeV)(MpL/lTeV) ~ 
10^^^ X (lTeV)(MpL/lTeV)^, where we have assumed that the decay is induced by dimen- 
sion four (three) operators for the first (third) expression. (The second one could appear 
accidentally.) The precise suppression needed depends, of course, on the details of the model. 
But it is clear that one needs a huge suppression factor for the decaying dark matter, if one 
would like to explain the PAMELA/ATIC data within the framework of particle physics 



16, a 



20l |. It is also clear that the existence of such a small number can not be explained 



in a low energy theory. 



III. THE MODEL: RADIATIVE SEE-SAW AND DARK MATTER CANDIDATES 

Here we would like to supersymmetrize the model of [4]. (An first attempt has been 



made in |27|.) We assume the R parity invariance as usual. So, we have R x Z2 discrete 
symmetry at low energy. The matter content of the model with their quantum number is 
given in Table I. L, H^,H'^ and rj'^,rj'^ stand for SU{2)l doublets supermultiplets of the 



Since 6 > 0, C + TV should appear instead of C for a negative C. 
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TABLE I: The matter content and the quantum number. Z2L is a discrete lepton number. Z2 is 
a subgroup of Z4, which is assumed to be anomalous and spontaneously broken by VEV of S and 
a down to Z2. 



leptons, the MSSM Higgses and the inert Higgses, respectively. Similarly, SU{2)l singlet 
supermultiplets of the charged leptons and right-handed neutrinos are denoted by and 
N^. (f) is an additional neutral Higgs supermultiplet which is needed to generate neutrino 
masses radiatively. E and a are also additional neutral Higgs supermultiplets which are 
needed to derive the superpotential (jl]) from a Z4 invariant one. 

We first consider a x Z2 invariant superpotential below. Later on, using S and cr, we 
will describe a possibility to obtain it from a i? x Z4 invariant one: 

W = W4 + W2, (4) 

where 

= + F,^L,Ar|r/" + A„r/"//'^0 + Arfr/'^//> + ^j^i/"//^ (5) 

W2 = ^-^N^N^, +f^,r^y + \^,,<p\ (6) 

The Yukawa couplings of the charged leptons can be assumed to be diagonal without 
loss of generality. 

Soft-supersymmetry breaking terms are necessary to generate neutrino masses radiatively. 
For the relevant Higgs sector they are given by 

CsB = -ml^fj''^'' - mlafj''^'' - ml^^ - (Br^fi'^fi'' + h.c.) 

-(^5^02 ^ ^ [AuKfi^'H''^ + AAdV'^H''^ + h.c), (7) 

where the hatted field is the scalar component of the corresponding superfield. The B and 
A soft terms are responsible for the radiative generation of the neutrino masses. We assume 
that fj^^T]'^ and do not acquire vacuum expectation values. 

To calculate the one-loop neutrino mass matrix, we treat the B terms as insertions. 
Then we find a one-loop diagram with one insertion of i?^r)"r}'^, which mixes 17" and f}'^. 
Correspondingly, we define the approximate mass eigenstates r/^ (the neutral component of 
rf) as 

?7^\ ^ / cos^ sin^A ^g>j 
?7q / I — sin cos 9 j \Vo I ^ 
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where 

tan2e = ml = i {ml + ml, ± [(m^ - ml)' + AmU'/'} , (9) 

m-uu m^ z 



with 



mid 



= + - ]^Ml cos 2/5 + ^A^t;^ cos^ /3, (10) 
= /i? + m^, + cos 2/3 + ^A>2 sin^ /?, (11) 
"^L = ^A„Arft;^cos/3sin/?, (12) 

and tan/3 = , v' = v\ ^ = {g^ + fl''^)'i^^/4. Neglecting higher order insertions 

we obtain the neutrino mass matrix at one loop: 

(M,),, = -^Y^,UiuMuUl^Y^^B^siYi2e[~cos'ei{m+,m+,Mu) 

+ sin2 6'/(m_,m_,Mfc) + cos26'/(m+, m_, M^)] , (13) 

where [/ is a unitary matrix defined by {U'^ Mis[U)ik = MkSik ^, and 



I{ma,mb,mc) 

"1 pl—x 

dx / dy[mlx + mly + ml{l — X — y)]^^ 



Jo 

^ mlml\n{ml/ml) + mlml\n{ml/ml) + mlml ln(mg/mg) 
(m2 - ml){ml - {ml - ml) 

As we can see from (jO]), (fT2l) and ( fT3l) . the neutrino masses are proportional to 5^ and A^A^ 
at the lowest order, because sin26' oc AuA^. So, the neutrino masses can be controlled by 
these parameters along with the Yukawa couplings Y^i, the masses of the inert Higgses and 
right-handed neutrinos. 



There are many candidates for the dark matter in this model [27|. The lightest combi- 
nation of each row in Table II could be a dark matter. But there can exist only three types 
of dark matter including the left-handed neutrinos depending on which discrete symmetry 
guarantees their stability. We assume that the first right-handed neutrino ip^i is the lightest 
one among ^/'at's and denote it by ipN (its mass is denoted by mNDAi)- So, ipN and the 
lightest neutralino x ( its mass is denoted by m^DM) are the dark matter candidates. Both 
have an odd R parity, so that one of them can be the stable dark matter, while the other 
one is the decaying dark matter, if Z2 is broken. 



Following 39|, |40t], we have computed the thermally averaged cross section for the anni- 



hilation of two V^at's and that of two x's by expanding the corresponding relativistic cross 



Ml < M2.3 is assumed, and we denote Mi by mpfDM later on. 
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TABLE II: The dark matter candidates. The (+, +, — ) candidates are dropped, because they are 
the left-handed neutrinos. 

section a in powers of their relative velocity, and we have then computed the relic densities 
^NDM and fl^DM- We have assumed that the SM particles are the only ones which are 
lighter than ipN and x, so that we have used the SM degrees of freedom at the decoupling, 
i.e. (7* = 106.75. We have found that, for the given interval of the dark matter masses, i.e., 
1 TeV ^ rrij^DM ^ 3 TeV and 0.2 TeV ^ rrixDM ~ 0.5 TeV, there is an enough parameter 
space in which {Qndm + ^xDM)h^ — 0.11 is satisfied. In the next section we let ipN decay 
into Xi while emitting high energy positrons. As it is clear from the superpotential ([5]), iIjn 
can not decay into the quarks, because ?7's do not couple to the quarks. 

IV. DECAYING RIGHT-HANDED NEUTRINO DARK MATTER AND 
PAMELA/ATIC DATA 

As long as the discrete symmetry i? x Z2 is unbroken, there are two CDM particles in 
the present model. One finds that the R x Z2[SU{?>)cf' and R x Z2[SU{2)l\^ anomalies are 
canceled with the matter content given Table 1 ^. Our assumption is that Z4 is anomalous 
and spontaneously broken to its subgroup Z2. Note that forbids W2 in ([6]) while W4 is 
allowed. Therefore, we have to produce it from an additional sector. This situation can be 
realized as follows. Consider the Z4 invariant superpotential including the SM singlet S and 
a given in Table 1: 

+ ( ^^ATf at; + A,r/V + ^A,0^) S. (15) 

The superpotential f|T5|) serves for S and a to develop VEVs, and consequently, Z4 is spon- 
taneously broken to Z2, producing effectively the superpotential IQ. The true stable dark 



^ We do not consider the R x Z2[U{1)yY ^'^d mixed gravitational anomalies, because they do not give us 
useful informations. 



7 



matter is the lightest one which has an odd parity of R. In the following discussion we 
assume that ^/'at is heavier than x- Since the ATIC data are indicating that the mass of the 
decaying dark matter particle is preferably heavier than 0(1) TeV, all the superpartners 
should be heavier than 0(1) TeV if m^DM > ^ndm- It is, therefore, more welcome for 
ipN to be the decaying dark matter, because a heavy ipN means a heavy r] Higgs, which is 
desirable to suppress FCNC processes such as /i ^ e + 7. 

As one can find, Z4 is anomalous: AQ = 1 mod N/2{= 2) ^ . Consequently, the suppres- 
sion coefficient 6 of ([3]) can take values 

Att^C , 2 C (mod 4) 

b = -TT^ = 47r2 X } — ^, (16) 

AQ 1 (mod 2) ' ^ ^ 

where C is the charge of exp{—bS). We assume that the non-perturbative effect can generate 
R invariant, but Z4 violating operators. At d = 3 there is only one operator rj^L which is 
even under R, and has the Z4 charge one. So we focus on t]^L: 

Wk = ^luV'^Li with/XM = P*MpLe-'<^\ (17) 

where pi are dimensionless couplings. Since (Fs/f) ~ ""^3/2 and {ip) ~ 0(1), the superpo- 
tential Wh induces a soft-supersymmetry breaking term 

Cb = Buff Li with Bu = wpjMpLmg/ae"'' (18) 

at the Planck scale, where is w a dimensionless constant. Since the Z4 charge of is 1, the 
charge of exp{—bS) has to be —1 mod 4, and then 

6 = 47r^ X (■ ■ ■ 7/3, 11/5, 11/7, 7/5, I---), 

which could give a huge suppression factor. 

With this observation we proceed with our discussion. The tree diagrams contributing 
to the tpN decay are shown in Fig. [1], where we have assumed that x is the pure bino. We 
do not take into account the tree diagrams which exist due to the mixing of ipr^u and tpe^, 
because these diagrams are suppressed by mf/rrir^, where m/ is the lepton mass. So, in 
the lowest order approximation only dimension two operators in the B soft-breaking sector 
exist. At the lowest order, ip^ can decay only into the leptons along with a x- parity 
violating operator LH^ allows the decay into the quarks, too.) The differential decay width 
is given by 



dTe+ "^NDM 2 f 1 



dE 7687r3 V 1 - 2x 



xM 1 - 



2 



Ai{l -2x- z') + 2Ai(l - x)(l + 2- 



z 



2 



2x' 



+6^2^ + 6^3(1 - 2x) , (19) 



For the Green-Schwarz cancellation to work, the Z/^SU (?>)c^ anomaly has to be matched to Z^\S\J (I^l^ 
anomaly. To realize this we introduce, for instance, a pair of 3 and 3 of SU (3)c with the charge one. 
Their mass can be obtained from < S > 3 x 3. 
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FIG. 1: The diagrams contributing to the ipN decay. We assume that x is the pure bino. The emitted 
positrons are all right-handed. Similar diagrams exist because of the mixing of 77" with E'^. The amplitude 
is suppressed by e~''m//m,,, so that we do not consider them. 



and the total decay width is 

Te^T = r^'oM = (1 - + ^3)(1 - 7.^ - 7z^ + z 



122887r3 

+4^2^(1 + 10^2 + /)] + 2Az'[-{Ai + A:i)z + 2^2(1 + z^) \nz\ }, (20) 



where 



^ = ^^^<1, x = ^—<{l-z')/2, (21) 

rriNDM ^NDM 

A, = 2g'' Yl l^/i^^ r ^3 = (22) 

^2 = -g" J2 KYnB*mB,) + h.c] (23) 

where g' (~ 0.345) is the U{1)y gauge couphng constant (the bino is assumed to be x)- 3 
runs over the negatively charged leptons, and i stands for a positively charged lepton in (122|) 
and (!23|) . We have assumed that all the scalar partners of the left-handed superpartners 
Il have the same mass m^. The positron can come from the decay of the anti-muons and 
anti-taus. In the following calculations, however, we assume that the energy spectrum of 
the positron coming from the anti-muon and tau does not differ very much from that of the 
direct production of the positron. So, we also sum over i = e+,/i"'",r+ in fl22l) and fl23|) to 
obtain dTg+fdE. At this order, all the emitted positrons are right-handed, as one can see 
from Fig. [H 

Before we calculate the positron spectrum, we briefly consider the suppression factor we 
need for our case. Assuming that Yij ~ 1 and that all the pi in 5's in (fT8l) are of the same 
size, we obtain 



/ TeV ^^ f mlrhl \ f tjindm ^ ^ 
tndm ~ I 

\mNDM 



) ( T'"^' 1 { T^^^l (P.^)-' (lO-^e^^) X 10^^ sec. (24) 
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FIG. 2: (e+/(e+ + e ) versus the positron energy E. The blue hues are the predictions of the model, 
where we have used: z = l/5(dashed), 1/6 (dot-dashed), l/5(dotted), TNDM{0-^^/t^NDMh^) = 1.4(dashed), 
2.0(dot-dashed), 3.0(dotted) x lO'^^ sec, niNDM = 2.0(dashed), 1.5(dot-dashed), l.O(dotted) TeV. The 
red points are the PAMELA data 23], where the predictions are written over the figure 4 of the PAMELA 



paper [23|. The solid line is the background published in 22], and it agrees with the one calculated from 
(P7)) - (|29l) without the primary source of the positron. 



So, we have a right order of t^dm with b = 47r^(7/3) which gives a suppression factor of 
10-8° (see (fT6|)). 

Now we come to compute the positron spectrum: 



dE'Ge+{E,E' 



dE' 



(25) 



where E'^ax = {rn%jj^.^ - m\j^^j) /2mNDKh dVe+{E)/dE = {t^dm) ^dn^{E)/dE, and we 
vary t^j:,m freely to fit the data. The positron Green's function Ge+ of |41| can be approxi- 
mately written as 13 | 



Ge+iE, E') 



0.11 



10 



16 



exp[a + b{E^-^ - E"-^)] cm^^ s 



1/5- 1\ 



E^ 



(26) 



where a^,6 depend on the diffusion model |13l . |42| . |43[ |. Here we use those of the MED 
model |43|]: a = —1.0203,6 = —1.4493,5 = 0.70, and we have assumed that except for the 
normalization factor Qj^i^Mh"^ /O.ll the decaying dark matter t/'^ has the same density profile 
in our galaxy as the NFW profile 4^. The background differential flux for each species are 
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FIG. 3: The differential energy spectrum scaled by E^. The solid, dashed, dot-dashed and dotted blue 
lines are the predictions of the model. The parameter values used here are the same as for Fig. [51 The solid 
blue hue is calculated with z ~ 1/10, Tjv£)A/(0.11/w7V£)M/i^) = 0.94 x 10^^ sec and itindm = 3 TeV. The 



predictions are written over the figure 3 of the ATIC paper |23| . where the PPB-BETS data j46| are also 
plotted. The black dashed line is the background presented in 23|. The normalization factor = 0.76 in 
((27| - (|29p is so chosen that the background computed from ((27|) - (l29l) agrees with the black dashed line at 
low energy. 
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^prim.bkg^^^ = N^O.IGE-'-' [1 + + 3.2^2-15] , (27) 

^sec.bkg^^^ = Ar<^0.7E°-^ [1 + 110^1-^ + 600^2-9 + 580^^-2] , (28) 

^sec.bkg^^^ = Ar^4.5E°-^ [1 + 650^2-3 + 1500E^-']"' (29) 

in the units in [GeV cm^ s sr]"^, where the energy E is in the units in GeV, and A^^ is 
a normahzation factor which we fix to be 0.76 from the ATIC data at low energies. The 
primary positron differential flux <|)|^™' is (c/47r)/e+, where fe+ is given in (l25l) . Using 
fl25|) - fl29|) . we then calculate appropriate quantities for PAMELA and ATIC: 

+ ^prim. ^sec.bkg 

— = : ^ for PAMELA, (30) 

g+ _j_ g- ^prim. , ^sec.bkg ^prim.bkg ^scc.bkg ' V / 

e+ e+ e~ 

= ^3^^prim. ^ ^sec.bkg ^ ^prim.bkg ^ ^sec.bkg^ j^rj.^^ ^3^^ 

dE e e e e 
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The results are shown in Figs. [2] and [3], where we have assumed Ai = A3 = A2 in (l22l) and 
f l23]) . The blue lines are the predictions of the model, and we have used: z = 1/5 (dashed), 
1/6 (dot-dashed), l/5(dotted), TNDAi{0.ll/ujNDMh'^) = 1.44(dashed), 2.0(dot-dashed), 
3.0(dotted) x 10~^^ sec, tundm = 2.0(dashed), 1.5(dot-dashed), l.O(dotted) TeV, where 



'221 



z is defined in fl2Tl) . The predictions are written over the figure 4 of the PAMELA paper 
and the figure 3 of the ATIC paper |23|. We see from Figs. [3] that ttindm should be heavier 
than 0(1) TeV in this model, too. 



V. CONCLUSION 

We have studied a dark matter model, in which one decaying and one stable dark matter 
particles coexist. We have assumed that one of the discrete symmetries ensuring the stability 
of the dark matter particles, when imbedded into a larger group, is anomalous, and the 
heavier dark matter can decay non-perturbatively. The huge suppression factor for the 
decay of dark matter to be needed can be obtained in this way. The concrete model we 
have considered is a supersymmetric extension of the Ma's inert Higgs model, so that the 
decaying dark matter (the lightest right-handed neutrino) can decay only into leptons along 
with the stable dark matter (LSP). We have shown that this scenario can explain the data of 
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23| . It is clear that if the recent and future data coming from the cosmic ray observations 
are intimately related to the nature of dark matter, its explanation may open the window 
to new physics beyond the SM. The radiative dark matter decay and high energy neutrino 
productions will be our next projects. 
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